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TECHNICAL NOTE 2771 


THERMAL BUCKLING OF PLATES 

By Myron L. Gossard, Paul Seide, and 
William M. Roberts 


SUMMARY 


An approximate method, based on large-deflection plate theory, for 
calculating the deflections of flat or initially imperfect plates sub- 
ject to thermal buckling is outlined. The method is used to determine 
the deflections of a simply supported panel subjected to a tentlike 
temperature distribution over the plate surface. Experimental results 
for a particular panel are in good agreement with the theoretical 
results for the range of temperatures and deflections considered in 
the test. 


INTRODUCTION 


In supersonic aircraft, deflections of plate elements out of the 
plane of the plate may be caused by aerodynamic heating. Deflections 
may occur without thermal stresses appearing if the plate is unre- 
strained and if the temperature distribution is linear throughout the 
volume of the plate. If the plate is restrained or if the temperature 
distribution is nonlinear, however, thermal stresses are induced. 
Deflections of the plate occur at the beginning of heating if the tem- 
perature varies throughout the plate thickness but do not appear until 
a critical temperature is reached if the temperature is constant 
throughout the plate thickness. This last type of distortion is an 
example of buckling of a flat plate by middle-surface forces that vary 
throughout the plate and is the subject of the analysis of the present 
paper. 

Since buckling of aerodynamic surfaces may have an adverse effect 
on aircraft performance, the temperature distributions for which buck- 
ling will occur should be known. In many cases it may not be feasible 
to design a structure which will not buckle and for these cases the 
magnitude of the distortions should be known in order that their effect 
on aircraft performance may be estimated. 

For initially flat plates, solution of the Von Karman large- 
deflection equations for plate buckling, modified to take into account 
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the effect of thermal stresses, is required for the determination of 
the buckle magnitude. A procedure for the approximate solution of these 
equations is outlined in the present paper and is illustrated by the 
determination of the buckle magnitude of a simply supported plate that 
is subjected to a tentlike temperature distribution. This temperature 
distribution was chosen because it was easily obtained experimentally. 

An approximate method for readily extending the results for initially 
flat plates to plates with initial imperfections is also presented. 

The effects of plasticity and of variations in material properties due 
to temperature are excluded from the analysis. 

The over-all validity of the assumptions made in the present analy- 
sis is tested by experimentally determining the variation with tempera- 
ture of the deflections of an initially imperfect plate subjected to a 
tentlike temperature distribution and comparing these experimental 
results with the theory. 


SYMBOLS 


a 

®mn 

b 

c l* c 2> * * * 
D 


D* 

E 

E* 


half-plate length in x-direction 

coefficients in series expansion for plate deflection 
half -plate width in y- direct ion 
coefficients in stress function F^ 


plate flexural stiffness, 


Et 3 


12(1 - ji2) 

transformed plate flexural stiffness. 
Young's modulus of plate material 


E*t* 3 


12(1 - \i 2 ) 


transformed Young's modulus of plate material. Ell 


W ir 


W_ 


F stress function defining stress distribution in plate, 

F = Fq + Fj_ 

Fq stress function defining thermal stresses in unbuckled plate 
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m ,n,p,q 


b 2 Eoffi 0 't 


b^EalOcr^ 

it 2 !) 


stress function defining stresses due to stretching of 
plate middle surface during bending 

plate thickness 

transformed plate thickness , t 
integers 

temperature distribution in plate 

temperature differential, difference between center and 
edge temperatures in a tentlike temperature distribu- 
tion (see fig. 2) 

critical value of Tq 

uniform edge temperature (see fig. 2) 

potential energy of an initially flat buckled plate 

plate deflection 

initial deflection 

initial deflection at plate center 

center deflection 

plate deflection shape given by small-deflection theory 
coordinate axes (see fig. l) 

temperature-differential parameter 
] 

critical- temperature-differential parameter 

plate aspect ratio 
coefficient of thermal expansion 
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a* 


transformed coefficient of thermal expansion. 



?xy 

€ X ,€y 


0 X , CTy 


xy 


^o^Yo^xyo 


shear strain in plane of plate 

normal strains in plane of plate in x- and y-directions, 
respectively 

Poisson’s ratio of plate material 

normal stresses in plane of plate in x- and y-directions, 
respectively, positive for tension 

shear stress in plane of plate 

thermal stresses in unbuckled plate 


V 2 

V 1 * 


differential operator, 


dx 2 dy 2 


differential operator. 


= - 'r- 

bx^ bx^'byb by^ 


ANALYSIS 

Statement of Problem 


The studies of thermal buckling presented herein are made for the 
panel shown in figure 1. The panel is heated along the longitudinal 
center line by a uniform line source of heat and cooled along the edges 
by two uniform and equal line sinks of heat. This arrangement supplies 
temperatures in the plate which are constant through the thickness and 
distributed in a tentlike manner over the faces as shown in figure 2. 

All edges of the panel are restrained in a direction normal to the plane 
of the plate by simple rigid supports but are free to slide in the plane 
of the plate. 

The investigation is made in steps corresponding to the stages 
through which an initially flat plate passes as temperature is gradually 
increased. First the thermal stresses at temperatures below the critical 
are determined; next the critical temperature is found; and then the 
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behavior of the plate at temperatures above the critical is obtained. 
The effects of initial imperfections throughout the entire temperature 
range are considered in a concluding step. 

The procedure is presented in the succeeding sections together 
with pertinent results of the mathematical study, the details of which 
are given in appendixes. Specific numerical results are given for a 
panel with aspect ratio a/b equal to 1.57 which is the aspect ratio 
of the panel described in the section entitled "Experimental Results." 


Thermal Stresses in the Unbuckled Plate 


Details of the calculation and experimental verification of the 
thermal stresses in the unbuckled panel are given in reference 1. The 
calculations of reference 1 employ the first-order approximation that 
on any cross section normal to the x-axis the stress a x is distributed 


Fq can then be 


as shown in figure 3. The thermal- stress function 

expressed as the product of a known function of y and an arbitrary 
function of x. The principle of minimum complementary energy is then 
used to determine an ordinary linear fourth-order differential equation 
for the function of x. The resulting approximate expression for Fq 


is 


0 - 15 b 2 BctT 0 (l - 3 ^1 + 2 4 )(b x ei 

H 


inh R^ ^ sin Rp ^ + 


Bp cosh R]_ — cos Rp 


(b < y < 0) (1) 


where Bp, R^, and Rp are defined in appendix A. 

The stress distribution in the plate is given by 


% = 


fc 2 


- EccTq^ - sinh Rj_ ^ sin Rp & + Bp cosh R]_ ^ cos Rp ^ + 1^ 


(2a) 
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0y O = Sx? 






sinh Rj_ ^ sin R 2 ^ + 


D 2 cosh R^ — cos R2 

Si 


i) 


K 'Jo ' ax 3y 


(2b) 



( l 


sinh R x ^ cos % f + % cosh B 1 | sin R 2 “) 


(b < y < 0) (2c) 


where Di, I) 2 > D3, & nd are defined in appendix A. The stresses in 
the region -b ^ y 5 ? 0 are identical with those given by equations (2). 

The stresses thus are a function of the temperature differential Tq 
for the tentlike temperature distribution and are independent of the 
edge temperature Tq_. 


Critical Temperature Differential 

As the temperature differential Tq of the heated plate increases , 
a value To cr is reached at which the plate buckles under the action of 

the induced thermal stresses. If only small deflections of the buckled 
plate are considered, the assumptions can be made that the middle sur- 
face of the plate does not stretch and hence that the stress distribu- 
tion in the plate does not change after the onset of buckling. The • 
stress distribution then is given by equations (l) and ( 2 ). The deflec- 
tion of the buckled plate is governed by the differential equation 
(ref. 2) 
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D7^w = t 


V x o + 


d 2 w ^ c^w \ 

% Sy 2 x yo a x a 7 J 


(3) 


and the critical temperature differential may he found ty methods 
appropriate for the investigation of the stability of flat plates with 
internally varying stresses. 

One such method is the Rayleigh-Ritz energy method which employs 
the principle of minimum potential energy (see ref. 2). For the present 
problem a buckle pattern symmetrical about the center of the plate is 
chosen as 


w 


CO 00 

51 51 Smn cos 

m=l,3,5 11=1,3,5 


mroc 

2a~ 


cos 


nity 

2b 




which, together with equations (2), is substituted into the potential- 
energy expression (ref. 2) to yield 




(5) 


Equation (5) is then minimized with respect to the unknown coeffi- 
cients amn. This procedure leads to a set of simultaneous equations 
which constitute a characteristic-value problem, the solutions of which 
give sets of relative values of the coefficients a^ and associated 
values of the critical temperature. The simultaneous equations are 


b 2 EaT n t 
u cr 


KpqSpq + ^ s* Kpqmn^n - 0 

m=l,3,5 n=l,3,5 


jt^D 


(p — 1, 3, 5, • • • > 
q = l, 3, 5, . . .) (6) 

Equations for the coefficients K pq and K pqmn are given in appendix A. 
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One method of solving equations (6) is a matrix iteration process 
which is described in reference 3- The lowest critical temperature 
coefficient thus obtained for a panel with aspect ratio of 1.57 is 




5-39 


(7) 


Since E appears in both the numerator and denominator (in D) of the 
critical- temperature-differential parameter^ the critical temperature 
differential is independent of Young’s modulus E. 

The approximate buckle pattern associated with this critical tem- 
perature is given by the equation 

“ = § 003 ft * °- 0365 C0B H cos lr + 

0.1360 cos ^ cos ^ + 0.0042 cos C os — ) (8) 

2a 2b 2a . 2b / 

in which w c is the deflection at the center of the plate and ^ = 1.57* 

A pictorial representation of equation (8) for one quarter of the plate 
is shown in figure 4. 


Post-Buckling Behavior 


In order that the post-buckling behavior of the heated plate may 
be obtained, stretching of the plate middle surface due to bending must 
be taken into account. Because of this stretching, the stresses in the 
plate change as the plate deflects and are determined by the equations 


V^F = -EaV^T + E 



8 2 w 8 2 w 

8x2 £jy2_ 


(9) 


°x = 


and 


8 2 f 

dy 2 


(10a) 



NACA TN 2771 


9 


CT y = 


d% 

dx 2 


*y 


dx dy 


(10b) 


(10c) 


in conjunction with the condition that the boundaries are stress free. 
The stress function is also related to the deflections by the equation 
of equilibrium 


= J&L ^3 * _ 2 a2p s2w N 

Vdy 2 dx 2 dx 2 dy 2 dx dy dx dyy 


( 11 ) 


so that equations (9) and (ll) must be solved simultaneously. 

Equations (9) and (ll) are the Von Karrnan equations (ref. 2) for 
large deflections of a plate , modified for the effects of thermal 
expansion (see appendix B). Exact solutions of these equations are in 
general difficult to obtain^ and approximate methods of solution must 
be used. A procedure vising the Galerkin method (ref. 4) is used in 
the present paper. 

The stress function F is obtained as the sum of two parts 

F = F 0 + F x (12) 

The stress function Fq is the thermal-stress function for the unbuckled 
plate which satisfies the equation 

V^FQ^-EctV^T (13) 

and the boundary conditions on stresses. This solution is given by 
equation (l). The stress function Fq is the solution (satisfying 
the boundary conditions on stresses) of the equation 




d 2 w d 2 w 
dx 2 dy 2 


(14) 


in which w is the buckle pattern given by equation (4). An approximate 



10 


NACA W 2771 


solution for Fq is obtained by the Galerkin method as set forth in 
reference 4. The stress function F]_ is chosen as the sum of a series 

of functions for which the boundary stresses vanish, suggested in 
reference 5 as 


F|_ = (x2 - a2)^(y2 - b2) 2 (c-^ + Cgx2 + c^y2 + . . . ) 


(15) 


The coefficients c lf c 2 , c^, . . . are found in terms of the coeffi- 
cients a^ by the equations 



(i = 1, 2, 3, . . .) ( 16 ) 

The resulting stress function F is substituted into equation (ll). 
The Galerkin method may again be used to determine the values of the 
coefficients a^ of the deflection function w. A set of simultaneous 
equations is obtained 



cos 


mjtx 

2a 


cos 


nity 

2b 



&2p c^w b 2 w 

Sy^ cix 2 cbc^ ciy 2 


2 


a% a 2 w \ 

dx Sy ax ay/ 


dx dy = 0 


(m - 1, 3) 5* • 
n = 1, 3, 5. . 


.) 


(17) 


which can be solved for values of the coefficients a mn . 

Equations (17) are nonlinear and their solution becomes difficult 
if many terms are retained in the deflection function. Experience has 
shown that very good results may be obtained if the shape of the 
deflected surface of the plate for large deflections is taken as the 
one existing at the critical temperature - that is, only the coeffi- 
cient ajj is left arbitrary and the ratios a^/a^ are assumed to be 

those given by the small-deflection solution previously described. 

The Galerkin equation from which the coefficient a^ can be determined 
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then becomes 



£ 2 F d^Vq S^Wg \ , 

— - — 2. + 2 - 2 -£ 2_ dx dy = 0 (18) 

dy 2 dx 2 dx 2 By 2 dx dy dx dy/ 


where 


Wg 


£ $: sg 

m=l,3,5 n=l,3*5 1 


mrac 

2a~ 


cos 


niry 

2 b - 


(19) 


and the values of a^/a-) -[ obtained from the small-deflection solution 

are also substituted into the stress function F. This procedure yields 
a relationship between the coefficient app and the temperature dif- 
ferential T 0 . For a plate of aspect ratio a/b equal to 1.57* the 
relationship found in this manner is- shown in appendix C to be 


b 2 EoT 0 t 

it%) 


5.39 + 1 . 12(1 - 



( 20 ) 


where w c is the plate-center deflection and is equal to 1 . 1767 a -) 1 . 

It can be seen that the plate deflections are independent of Young *s 
modulus E, since E appears in both the numerator and denominator 
of the temperature -differential parameter. 


Effect of Initial Imperfections 

Since actual plates are not usually flat, initial imperfections 
should be taken into account in the previously developed analysis. 

The analysis of the thermal buckling of initially imperfect plates 
involves the solution of the Von Kerman large-deflection equations 
for initially imperfect plates modified for the effects of nonuniform 
temperature distributions (appendix B): 

S 2 Wj_ \ 2 d^j^ d^.^ 

dx dy/ dx 2 dy 2 


V^F = -EaV 2 ! + E 


f d g w \ 
\dx dy/ 


d^ d^w 
dx 2 dy 2 




( 21 ) 
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D / \ S2w _ S2 f b 2 v 

£V^(w - Wi) + 2 

t by 2 bx 2 bx 2 by 2 dx &y ox by 


( 22 ) 


The solution of these equations could he effected in much the same 
manner as was done for the equations for flat plates. The outlined 
process, however, is tedious, and an approximate method for assessing 
the effect of initial imperfections would he more advantageous. Such 
a method is developed in appendix D. 


In order to analyze an initially imperfect plate, a flat plate 
having the same aspect ratio is first analyzed, the quantities E, a, 
and t being left arbitrary. Then everywhere in the resulting expres- 


w 


. 2> 


sions, E is replaced by E[l 
t 


w c 2 


, a by it, and t by 

1 1 y ic 2 


W r 


This relatively simple procedure yields the stresses and 


i » c 

deflections of an initially imperfect plate. 


For the problem of thermal buckling of a flat, simply supported 
plate of aspect ratio 1.57* subjected to a tentlike temperature dis- 
tribution, the relationship between the temperature differential Tq 
and the center deflection of the plate w c has been found to be 


b 2 EaT ( - ) t 

Tt 2 D 


5.39 + 1.12(1 - 



(23) 


When a and t (note that t also appears in D) are replaced by 
and ■ - - (no substitution need be made for E since 


v*t 2 
w ic 




w c ^ ’ ■'c 

the relationship is independent of E) the resulting equation may be 
written as 


b 2 EaTot 

ji^D 



+ 1 . 12(1 - 


V?) 


Wc 2 - Vi c 2 
t 2 


m 
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which, is the relationship between the temperature differential and the 
plate-center deflection for an initially imperfect plate. 

This method of analysis may be expected to be fairly accurate for 
plates for which the initial deflected shape is similar to the shape 
of the first buckling mode of the corresponding flat plate. If the 
first buckling mode does not predominate in the initial deflected shape, 
no recourse appears possible except to solve equations ( 21 ) and ( 22 ) 
for the large deflections of an initially deflected plate. 


EXPERIMENTAL RESULTS 


In order to check the applicability of the theoretical analysis to 
actual plates, the deflections of a simply supported rectangular plate 
subjected to a tent like temperature distribution were obtained experi- 
mentally and were compared with the theoretical results. The panel 
was tested under conditions which gave virtually complete boundary 
freedom in the plane of the plate and simple support. 

A plate (see fig. 5) having an over-all length of 3 6 inches, an 
over-all width of 2k inches, and a thickness of 0.25 inch was used; 
this plate is similar to the panel of reference 1. The length and 
width between simple supports were 35-25 inches and 22.50 inches, 
respectively, which correspond to the plate aspect ratio of 1.57 used 
in the numerical calculations previously reported. The coefficient of 
thermal expansion a and Poisson’s ratio p for the 75S-T6 aluminum- 

alloy plate material were 0.127 X 10“^ ■ —•y e — and O. 33 , respectively, 

in . - F ■ 

for the range of temperatures of the test. The plate was initially 
imperfect and had a center deflection of 0 . 0^5 inch. 

The arrangement and operation of the heat source and sinks are 
the same as for the panel of reference 1. Deflections of the plate 
along the longitudinal and transverse center lines of the sheet were 
measured by dial gages. These deflections are the difference between 
the total deflections and the initial deflections of the plate. 

Figure 6 shows the experimental and theoretical plate-center 
deflections w c - w-j_ c plotted as a function of the temperature dif- 
ferential Tq. The theoretical curve is given by equation ( 2k ) with 

P, a, t, b, and w^_ replaced by their respective values. With 

c 

a = 0. 127 X lO"^ in ’ h 
in . -°F 
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t = 0.25 in. 


b = 11.25 in. 


w i c = 0-0^5 in. 

and 

H = 0.33 

equation (24) becomes 

T 0 = 194.l(l - ~^) + 573. 6 (w c 2 - (0.045)5 (25) 

where Tq is in degrees Fahrenheit and w c is in inches. Figure 6 
shows that good agreement exists between theory and experiment for the 
range of temperatures and deflections considered in the test. 

Also shown in figure 6 is the theoretical variation of center 
deflection with temperature for an initially flat plate (wi c = 0). 

The additional deflection caused by initial imperfections of the plate 
can be seen to be appreciable. 

The deflections measured along the longitudinal and transverse 
center lines of the plate for various values of the temperature dif- 
ferential Tq are shown in figure 7. The theoretical deflections 
calculated for the same values of Tq are also shown in the figure. 
Although the measured and calculated initial deflection shapes (Tq = 0° F) 
are not in very good agreement, the measured and calculated deflection 
shapes tend to come into closer agreement as the temperature differential 
becomes larger and thus tend to fulfill the assumption of the theoretical 
analysis that the plate deflections have the same shape as the shape of 
the first mode of buckling of the corresponding flat plate. 


CONCLUDING REMARKS 


Good agreement exists between experimentally determined center 
deflections of an initially imperfect panel subjected to a tentlike 
temperature distribution and the center deflections calculated from 
an approximate solution of the large-deflection equations for an 
initially imperfect plate. Other cases of thermal bending or buckling 
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of plates can very likely "be calculated with good accuracy by the 
methods developed in the present paper if the temperature distribution 
and initial imperfections of the plate are known. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , May 12, 1952 
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APPENDIX A 


EQUATIONS FOR CRITICAL TEMPERATURE AND BUCKLE PATTERN 


The equations for the critical temperature and buckle pattern for 
symmetrical buckling of the panel subjected to the particular tempera- 
ture distribution considered herein are as follows (eqs. (6)): 


b 2 EaT r 


'cr 


npq^q + 2 — 2 — ^qmn^n " 0 

t m=l,3,5 n=l,3,5 


2 d 


(p = 1, 3, 5, 
q = 1, 3, 5, 


. . ) (Al) 


The coefficients Kp q and Kp gmn are defined as follows: 

,2 


Ki 


■pq 


^ I * 


(A2) 


%qmn - P |roAnq^®lAnp + ®2®mp + ^mp) + ^^nq^S^P + ^^mp^J 
q |^nq( D l!mp + ^Jmp) + mBqn^Gpm + D^Hp^yj 


The constants in equation (A3) that apply also in the formulas for 
direct thermal stress (eqs. (l) and (2)) are given by 


(A3) 


k-i sinh R-i cos Rp - k; 
B X = — - - 


/■» 


(A4a) 


k^ cosh Rq sin Rg - kg sinh R^_ cos Rg 
Bg = - ■■■■ ■ \ AAb ) 

kj_ sin Rg cos Rg + kg sinh R]_ cosh R-j_ 
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D 1 = B l( k l 2 “ *z) ~ 232*1*2 
D 2 = B 2 (k x 2 - + 2B 1 k 1 k 2 

D 3 “ B l k 2 + ^*1 
= B]_k]_ - B 2 k 2 

where 


R 1 = k l t 
R 2 = *2 | 



The other coefficients are 



(A4c) 

(Akd) 

(A4e) 

(A4f) 




> (A5a) 


is even or zero 



18 


NACA TN 2771 




NACA TN 2771 


19 



sinh Rt — sin Ro — sin 
x a ^ a 


mitx . pitx 
7: — sin £ — 
2a 2a 


a 


m-p 


2 I — 

= (-1) KCi + C2 + + C^)sinh R^ cos R2 - (C^ + Cg + 

C7 + Cg)cosh R-^ sin Rgj 


(A5d) 




-r 


cosh Rt — cos Rp — sin SZ2£ sin EE a — 

A o Oo n 


2a 


2a a 


m-p 


(-1) 2 + C2 + Cg + Cj^)cosh R-j^ sin R 2 + (C^ + Cg + 

Cy + Cg)sinh R-^ cos R^l 


v- W 

II 

J* 

to 

->a 

sin 

) 

mitx 

2a 

, pitx 
sin i — 
2a 

— 0 

if 

m - 

_P 

+ 0 

x mp 

2 


F = — 

if 

|m - 

p 

= 0 

mp 2 

2 1 


> 


J mp 


sinh R-j f cos Ro | cos sin EE d £ 
a a 2a 2a a 


m-p 


= (-1) 2 [Tc, - C 2 + - C^)sinh R-j. cos R 2 - (C^ - Cg + 

C7 - Cg)cosh R^ sin R^J 


(A5e) 


(A5f) 


(A5g) 



20 


NACA TN 2771 


rr / , _ x . _ x mitx pitx , v 

= / cosh Rp— sin R 2 — cos - — cos £ — d£ 

^ J 0 x a 2 a 2 a 2 a a 


m-p 


(-1) 2 J^C-^ - C 2 + C 3 - C^)cosh Rp sin R 2 + (C^ 
Cj - Cg)sinh Rp cos R^j 


- c 6 


L mp 


-r 


. , -r. x . m x mrtx pitx , v 

smh Ri — sin Ro — cos cos — — d — 

A a a Oo Oo 


2 a 


2 a a 


(A5h) 


= (-1) 2 JT Cl + C 2 - - C^)sinh Rp cos R 2 - (C^ + Cg - 

Cy - Cg)cosh Rp sin Rgj 


(A5i) 


J mp - f 8 " cosh Rp — cos R 2 cos cos £52. d2 
V J 0 a a 2 a 2 a a 


^ +1 

= (-1) |7cp + C 2 - - C^Jcosh Rp sin R 2 + (C^ + Cg - 


C 7 - Cg)sinh Rp cos r£| 


(A5J) 



m + p 

R 2 + — 5 — * 



m + p 

+ 



where 
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If only the terms a-^, a^, a 31> an( ^ a 33 are retained in the 

deflection function (eq. (4)), equations (Al) may he written in matrix 
form, for a plate having an aspect ratio of 1.57> as 


15.73 

22.52 14.88 

-5.96 


c 

a ll 



a ll 

0.504 

1.426 0.871 

0.377 

< 

a 13 

> 

100 ; 

a l3 

1.35 

3.54 7-42 

7.69 


a 3l 


b 2 EoT 0 t 
u cr 

a 3l 

-0.0735 

0.208 1.043 

0.437 


a 33 

L. 


n 2 D 

_ a 33. 


(A6) 


The solution for the largest value of 10C> 

h 2 EoT 0 -t 


D 2 EaT 0 t 
u cr 




, and hence for the 


smallest value of. 




J cr 


f is obtained from matrix iteration (ref. 3) 


of equation (A6) . The critical-temperature parameter so obtained is 


b^Ealn t 

u cr 


75 " 


= 5-39 


(A7) 


The relative values of the four coefficients retained in the deflection 
function are found to be 


a ll 


r 

1 

J a 13 

> = « 

0.0365 

a 3l 


0. 1360 

a 33 


0.0042 


(A8) 


in which case the deflection function may be written as 


w 


= an (cos — cob + 0.0365 


itx 3ny - 3itx jty 

COS 7 r— COS + 0.1360 COS ^ COS + 

2a 2d 2a 2D 


0.0042 cos cos 

2a 2b J 


(A9) 
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The deflection at the center of the plate Is 1. 1767a]_]_. Let this 
quantity he denoted hy w c ; then 

w = — — — /cos — cos — + O.O 365 cos — cos + 
I.I767V 2a 2b 2a 2b 


O.I36O cos cos — + 0.004-2 cos cos (A10) 

2a 2b 2a 2b / 


which is the form of the deflection function presented as equation (8). 

The choice of the particular coefficients a-Q, a^, a^, and a^g 

depends on the fact that these are the most important coefficients in 
the series for the deflection w (eq. (4-)). The following table shows 
the convergence of the critical-temperature parameter as more terms are 
taken in the deflection function: 


Terms retained 

b 2 EctT 0 t 
u cr 

jt^D 

a ll 

6.35 

a ll> a 31 

5.65 

a ll> a 31' a 13 

5.40 

a ll> a 31> a 13> a 33 

5-39 


The retention of terms other than the four chosen has a negligible 
effect on the critical temperature and on the buckle pattern. 
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APPENDIX B 
/ / 

MODIFICATION OF VON KARMAN LARGE-DEFLECTION EQUATIONS 
FOR EFFECTS OF THERMAL EXPANSION 


The condition of compatibility of strains for a deflected plate 
with initial imperfections is represented by the following equation 
(see refs. 2 and 6): 


d 2 ex + d 2 e x 

Sy2 Sx2 



d 2 w d 2 w /"S^wi \ 2 S^wi S 2 wi 

Sx2 c)y2 \Sx by/ Sx2 Sy2 


(Bl) 


When heat is applied to a plate, strains e x , 6y, and are 

caused by thermal expansion of the material as well as by stresses that 
may arise from various sources. The strains are given by 


e x = |K - W y ) + ctT 
e y = |(ffy - H«r x ) + aT > 


7 


xy “ 


2(1 + n) 

E 



(B2) 


With the introduction of a stress function F defined by 


ff x " 


a y = 


'xy 


d 2 F 

Sy 2 

5 2 f 

dx 2 


> 


5% 

dx Sy 


(B3) 
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the strains are given by 






2(1 + ji) B 2 F 
7xy = - -SH- 


(*) 


Substitution of equations (B4) into equation (Bl) gives the compatibility 
condition in the form 


vV = -EkiV 3 ! + E 



52 y ^2w _ Y + 

Bx 2 ^y 2 V* W Bx 2 By 2 


(B5) 


This is the first of the Von Karman equations, modified for effects of 
thermal expansion. The stresses derived from the stress function F 
satisfy equilibrium in the plane of the plate because F, as defined by 
equations (B3), identically satisfies the equilibrium equations 


5t xy _ 
Bx By 


> 



bT. 


J2. 0 

Bx 


(B 6) 


The second Von Karman equation, the equation of equilibrium of forces 
normal to the plate middle surface, remains unchanged as 


W^(w - w^) 


'b 2 F B£v + 
py 2 Bx 2 


B 2 F B^ _ 2 
Bx 2 By 2 


B^F B^ 

Bx By Bx By, 


(B7) 



APPENDIX C 


ro 

CT\ 


/ f 

APPROXIMATE SOLUTION OF THE VON KAHMAN LARGE-DEFLECTION EQUATIONS 
Determination of the Stress Distribution In the Plate 


/ / 


The firBt Von Karman equation for flat plates relates the stress distribution in the plate 
to the plate deflection w as follows: 

v2 


vV = -EaV^T + E 


/jV\ d^w d^w 

W dy/ dx 2 by 2 


(Cl) 


The stress function F can he separated into two parts, the primary stress function Fq, which 
satisfies the equation 


vV n = -EcffT^T 


(C2) 


and is given by equation (l) for the present problem, and a stress function Fj_, which satisfies 
the equation 

7 \ 2 tA. a&T 

f V W j V w V w 

\5x byj dx 2 dy 2 


= E 


(C3) 


and the given condition of stress-free edges. For the buckle pattern given by equation (it-) , 
equation (C3) becomes 


1 * 1 -« £ £ £ 


m=l,3,5 n=l,3,5 p=l,3,5 q=X?3,5 


mit q* 
2a 


f it /pit nrt . mitx . pitx . nity . qity 

b 2b aln 2i- s “ ST sln 2^ ^ - 


nit qit mitx 

— -a- cos cos 

2a 2b 2a 


cos ss cos aaN 
2a 2b 2b J 


(C4) 
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Equation (C4) may "be solved approximately by the Galerkin method, A solution for Fq is 
taken as 


F 


1 “ 



c i f i 


(C5) 


Each function fq is a function of x and y that yields stresses which satisfy the given 
boundary conditions. The arbitrary coefficients cq are then obtained from the equations 


v g~ pa nb , pa Pb 

c q J J ^ f s V 4 fi dx dy = J J ^ f s G(x,y) dx dy 
where G(x,y) denotes the terms on the right side of equation (C4). 


(s - 1, 2, 3, ... ) (C6) 


For the present problem the boundary conditions to be satisfied are those of stress-free 
edges of the plate. Furthermore, the stress distribution is symmetrical about the x- and y-axes. 
A series for Fq that satisfies these conditions is suggested in reference 5 as 


Fq = (x 2 - a 2 ) 2 (y 2 - b 2 ) 2 (cq + C2X 2 + cgy 2 + . . .) 

In the present solution only the first three terms of the series are considered. 


(07) 


The substitution of equation (C7) into equations (C6) yields the following equations for 
the determination of the coefficients c^, c 2 , and c^* 


a 2 b 6 cq + (1 


( l + 4 A6= 2 + (4 + $) 


a^b^Co ■ - E 

3 32768 


r/i + is£ + ajS\ 

V 7 b 2 bV 

^13 2 + a 31 2 ) + (y + jF)( a 13 + a 3l) a li + + H+)( a 13 + a 3l) a 33 ’ y 


65§ 

45 


1863 

—r a i3 a 3i + 

TT * 


f(‘ 

162 


a ll 2 * a 33 2 1 + 


2 ) 


a ll a 33 + 


9 9 • 


(C8a) 
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e ♦ rr h $ * * ^ - - m • - iK + 

ftf - f)-3 2 * (SI - iSK * (f - ^ * (I + ^ - ^>-3 * 

(t ■ f + T - *&** - {f - if>-33 * (f - 5V + (f -| 


ro 

m 


189 3825Y /108 189 2id\ 

5" ' i^ - / 13 * 33 + VaT + X ' T/ 31 * 33 


4 - 


/ rt r>i \ 

^<jod; 


& ♦ sjw, ♦ i/i ♦ 1 3/Cfc ♦ 4aS * 4 W, , . . i§V i2 + 

\ XJ - hV - xx \ bV “ b* bV J 32768 |\35 n2/ XJ - 

/£72 . jA 2 + AMii . i*t) a 2 * (H* . i£l 2 + f36 . 90 + m . 3H2V + 

V 10 5 27*2/ 13 \ 315 t?) 31 \35 3 #) 33 \5 it 2 jt^ S J6 ) 11 13 

(y?&9 _ 60750^ , 


( 12 . 1701 1701^\ A86 6075\_ . 

ip^r- — j a n a 3i - - ^->^33 + + 


flSS. + 1SS. _ =“V„a„ + - H + i£ - 3tejV, n a„ + 

\35 *4 5 ^ ^ 2*&/ ^ ~ 


(C8c) 
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The assumption is now made that the buckled shape of the plate does 
not change as the deflections increase, in which case relative values 
of the coefficients a mn for large deflections are the same as those 
at the critical temperature. For the plate of aspect ratio a/b equal 
to 1. 57 , these relative values are given by equation (A 8 ). Only the 
coefficients a-j_p, a^g, a 3 i> a 33 are considered. The solutions 

of equations (C 8 ) are then 


C1 = -0.01337 


^ll 2 

*8 


C 2 = 0. 00685 


Ea ll 2 

*10 


> 


c 3 = 0.00389 


Ea ll g 

*10 


(C9) 


The stress function for stresses due to stretching of the middle 
surface of a plate with aspect ratio 1.57 is therefore given approxi- 
mately by 



0.00389 



(-1.57 ^ f ^ 1.57; 
-l<^l) (CIO) 


where w c is the plate-center deflection and is equal to 1 . 7 . 

The stress function F is given by the sun of equations (l) and (CIO). 


Determination of Relationship Between Temperature 
and Center Deflection 

With the stress distribution in the plate known ^as a function of 
the plate-center deflection w c , the second Von Karman large -deflect ion 



30 


NACA TN 2771 


equation for flat plates can be solved to determine the relationship 
between temperature and deflection for temperatures above the critical. 
Equation (B7) can be satisfied only on an average over the plate when 
the buckle pattern for temperatures above critical is assumed to have 
the same shape as at the critical temperature. Use of the Galerkin 
method (ref. 4) for this averaging process gives the following equation 
governing the relation between temperature and center deflection: 



a% a 2 F a% aV \ 

a y 2 a x 2 ax2 a y 2 Sx ^y Sx V 


0 (Cll) 


where 


V S 


£ X 


Qmn 

an 


mrtx 

cos — — cos 
2a 


my 

2b 


(C12) 


and the values of ajm/a^ are given by the small-deflection solution 
of appendix A. 


For the plate of aspect ratio 


a/b equal to 1.57, 


Wo = cos cos —■ + 0.0365 cos ^ cos 


Jty 


Itx 


2a 


2b 


2a 


Sr + °- 136 ° cos lr cos S + 


0. 0042 cos 


3«x 

2a 


cos 


3*y 

2b 


(C13) 


and the stress function F is given by the sum of equations (l) and 
(CIO). Equation (Cll) then yields the following relationship between 
the temperature Tq and the center deflection w c : 


b 2 EaT Q t 


5.39 + 1.12(1 - n 2 )^ 

t 2 


(Cl4) 
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The deflection at any other point of the plate with stress -free edges 
and of aspect ratio a/b equal to 1. 57 is given by 


w 

t 


0. 723 


-EoT 0 t 


- 5-39 


1.12(1 - p2) 


Wo 


(C15) 


where w g 


is given by equation (C13). 
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APPENDIX D 

CORRECTION FOR INITIAL IMPERFECTIONS 


The large -deflect ion equations for bending of an initially imper- 
fect plate subjected to thermal stresses are derived in appendix B as 


V^F = -EaV^T + E 



s vV - w t ) .. ** 

t 1 Sy 2 Sx 2 


5 2 w B 2 w j 

f \ 2 

B^.? 

j ' 


5x2 Sy 2 

^5x By y 

5x2 

By 2 

+ B 2 F 5 2 w 

o B 2 F 

B^ 


Bx 2 By 2 

Bx By Bx By 



(Dla) 


(Dlb) 


The assumption is made that the bending deflections of the plate are 
merely a magnification of the initial deflections w-jj that is, the 
deflections w have the same shape as the initial deflections w^. 
This as sumption may be written as 



w w c 


(D2) 


The substitution of equation (D2) into equations (Dl) yields 


vV 



c^w 5 2 w 
Sx 2 5y£_ 


(D3a) 



- 2 52f 

w c / 5y2 5x2 5x2 5y2 


(D3b) 


5x 5y 5x By 
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If the following substitutions are made: 



t* = 



equations (D3) become 


V^F 


-E*a*V%: + E* 



B 2 w B 2 w 
Bx 2 By 2 


m 


(D5a) 


E.r^ v = B^sfw + afw _ 2 &P B 2 w 

t* By 2 a x 2 Bx 2 By 2 dx By Bx By 


(D5b) 


Equations (D5) are identical with the large -deflect ion equations for 
buckling of a flat plate with Young’s modulus E*, coefficient of 
thermal expansion a*, and thickness t*. If the initial deflections w^ 
are assumed to satisfy the same homogeneous boundary conditions as 
would be satisfied by the deflections of an initially flat plate, the 
solution of equations (D5) is identical with the solution for the 
deflections of a flat plate having the same aspect ratio but with E 
replaced by E*, a by a*, and t by t*. 

As an example of this method of correction, consider the relation- 
ship between the temperature and the center deflection for the plate 
of aspect ratio 1.57 

b 2 Eo£ 0 t „ v 2 

— = 5.39 + 1.12(1 - n 2 )_|_ (d6) 

jt 2 !) t 2 
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To find the corresponding relationship for a plate with aspect ratio 1. 57 
hut having initial deflections of the same shape as the initial buckling 

deflections of the flat plate, a is replaced by a* 


and t (note that t also appears in D) by t* or 

(E need not be replaced by E* since the relationship is independent 
of E. ) Equation (D6) then becomes 



b 2 EaT 0 t 

*2d 



+ 1.12(1 - |i 2 )- 


v c 


2 _ 


w i. 


t 2 


(D7) 


In order that the validity of the foregoing method of analyzing 
plates with initial deflections be checked, the method was used to 
calculate curves of center deflection plotted against average edge 
compressive stress for the problem considered in reference 6, the 
bending under edge compressive stress of simply supported square 
plates with initial imperfections. These curves were in excellent 
agreement with the numerical results of reference 6 which are obtained 
from an approximate but accurate solution of the Von Karman large- 
deflection equations for initially imperfect plates. The agreement 
was found to exist for all cases in which the initial imperfection 
was a half-sine-wave deflection in both directions. 



NACA TN 2771 


35 


REFERENCES 


1. Heldenfels, Richard R. , and Roberts , William M. : Experimental and 

Theoretical Determination of Thermal Stresses in a Flat Plate. 

MCA TN 2769, 1952. 

2. Timoshenko, S. : Theory of Elastic Stability. McGraw-Hill Book Co., 

Inc., 1936 , pp. 305-323. 

3. Von Karman, Theodore, and Biot, Maurice A.: Mathematical Methods 

in Engineering. First ed. , McGraw-Hill Book Co. , Inc. , 1940 , 
pp. 196-203. 

4. Sokolnikoff, I. S. : Mathematical Theory of Elasticity. First ed. , 

McGraw-Hill Book Co. , Inc. , 1946, pp. 313-316. 

5. Timoshenko, S., and Goodier, J. N. : Theory of Elasticity. Second 

ed. , McGraw-Hill Book Co., Inc., 1951* p. 169 . 

6. Hu, Pai C., Lundquist, Eugene E. , and Batdorf, S. B. : Effect of 

Small Deviations From Flatness on Effective Width and Buckling 
of Plates in Compression. MCA TN 1124, 1946. 



36 


NACA TN 27T1 




Figure 2.- Tentl-ike temperature distribution. 
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Figure 4.- Small -deflect ion buckle pattern in one quadr an t of a plate of 

aspect ratio 1.57. 
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Figure 6.- Comparison of calculated and experimental deflections at 

plate center. 






